Exact solution of the Nonconsensus Opinion Model on the line 
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The nonconsensus opinion model (NCO) introduced recently by Shao et al., [Phys. Rev. Lett. 
103, 018701 (2009)] is solved exactly on the line. Although, as expected, the model exhibits no 
phase transition in one dimension, its study is interesting because of the possible connection with 
invasion percolation with trapping. The system evolves exponentially fast to the steady-state, 
rapidly developing long-range correlations: The average cluster size in the steady state scales as the 
square of the initial cluster size, of the (uncorrelated) initial state. We also discuss briefly the NCO 
model on Bethe lattices, arguing that its phase transition diagram is different than that of regular 
percolation. 
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I. INTRODUCTION 

A Nonconsensus Opinion model (NCO) has been intro- 
duced recently [l| where each node in a graph can be in 
one of two states (representing opinions). At each time 
step each node adopts the state agreeing with the ma- 
jority of the nodes in its neighborhood, consisting of the 
node itself a.nd its nearest neighbors. In the event of a tie, 
the node retains its original state 0- The NCO model 
is thus similar to the majority-voter model, but where 
self-opinion counts. This difference is sufficient to ensure 
survival of the minority opinion, in the steady state jH . 

In [l| a conjecture was made that the NCO model and 
invasion percolation with trapping (TIP) Q are in the 
same universality class, and this was well supported by 
convincing numerical evidence. As opposed to invasion 
percolation without trapping, which is known to be in 
the same universality class as regular percolation, TIP 
is much less well understood Q. The following exact 
analysis of the NCO model, even though limited to the 
simplest case of one dimension, is therefore interesting 
not only on its own right, but also because of the possible 
connection with TIP. In addition, we discuss the phase 
transition diagram of the NCO model on Bethe lattices, 
finding a rich three-phase diagram, reminiscent of regular 
percolation in nonamenable one-ended graphs. Models 
similar to the NCO have been considered in the context 
of social studies 1 41. 



II. THE NCO MODEL IN ID 

In the following we will specialize to the evolution of 
the NCO model on the infinite line. Each lattice site can 
be either empty (o) or occupied (•). The initial state is 
completely random, with a fraction = p oi occupied 
sites, and qo = I — p = q empty sites. At each time step 



all of the sites get updated simultaneously according to 
the NCO rule: the new state of each site agrees with the 
majority of the states of the site's neighborhood, consist- 
ing of the site itself and its two nearest neighbors. Thus, 
an occupied site remains occupied if at least one of its 
neighbors is occupied, but becomes empty in the next 
time step if both its neighbors are empty. An empty site 
remains empty if at least one neighbor is empty, and be- 
comes occupied if both neighbors are occupied. It is easy 
then to predict the density of occupied and empty sites 
after the first step: 

pi = p{l - q^) + qp" = 3/ - 2p3 , = 3q^ - 2q^ . (1) 

After this first step, the initially random state of the 
system develops correlations and one can no longer apply 
the same reasoning to predict the state of the system 
in subsequent steps. Indeed, ignoring correlations and 
iterating the mapping ([1]), one arrives at the erroneous 
conclusion that the final state of the system consists of 
either all empty or occupied sites (depending on whether 
p is smaller or larger than 1/2). 

It is clear from the evolution rules that a cluster of 
two or more occupied (or empty) sites remains stable 
forever. We denote occupied clusters •••••• with a solid 

rectangular box , and empty clusters o o • • • o with an 
empty box i — i . 

The initial configuration can be viewed as a collection 
of stable blocks, separated by alternating sequences of 
occupied-empty sites. There are only four possibilities: 



(a) I — I • o • o ■ 

(b) o • o • ' 

(c) I — I • o • o ■ 

(d) o • o • ■ 



(2n intervening sites) 
(2n intervening sites) 
(2n + 1 intervening sites) 
(2n + 1 intervening sites) 
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In all of these cases all of the sites in the alternating 
sequences flip state at each time step. On each step, the 
two sites at the edges merge with the bounding blocks. 
Thus, cases (a) and (b) cease to evolve after n steps, 
and the bounding blocks meet, having each grown by 
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n sites. For example, i — i • o • o^m — ^ i — i • o^m — >■ 

1=1^ . Cases (c) and (d) evolve for n + 1 steps, after 
which the bounding blocks merge, having subsumed all 
of the intervening sites, e.g., i — \ mom\ — r - 



III. THE STEADY STATE 

At the steady state the system crystalizes into alter- 
nating blocks, • • -c^^Mc^^M- • • , each consisting of at 
least two sites. We now compute several characteristics 
of the steady state. 

Let us first look at the density of kinks — the bound- 
aries between clusters. We consider kinks of type 
(kinks of type ^mc^ have the same density, by symme- 
try). 

The location of the kinks in the steady state can be 
predicted completely from the random initial state, as 
only case (a), above, develops into the type of kink in 
question. Thus the steady-state density of kinks. Kg, is: 
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Notice that just two adjacent occupied (empty) sites are 
sufficient to guarantee the existence of a ( i — i ) cluster, 
accounting for the overall p^g^ factor. That Ks is smaller 
than pq, the density of o»-kinks in the initial state, re- 
fiects the fact that the system coarsens over time. 

To obtain the steady-state density of occupied sites, ps, 
we observe that cases (a) and (b) keep the number of oc- 
cupied and empty sites constant over time. Case (d) con- 
verts n+1 empty sites into occupied sites, while case (c) 
does the opposite, emptying n+1 occupied sites. Putting 
all this together, we get 



Ps = p + p'*g(l + 2pq + ip'^q^ 



+ •••) 

-q^p{l + 2pq + ip^q^ + • • • ) 
pq{p'^ - q^) 



{l-pqf 



(3) 



Thus, unless p — 1/2, the majority species increases on 
expense of the minority species (Fig. [T]) . 

We now compute the probability distribution for 
clusters of n occupied sites in the steady state. In 
the initial state, the probability of an n-cluster is 
P{n) — q^p^, so the expected length is (ri,)o = 

j:ZinPin)/j:ZiP(n) = l/q. 

Fig. [2] lists all the possible scenarios leading to n- 
clusters in the steady state. Ultimately, a cluster of occu- 
pied sites needs be bounded by empty clusters on the left 
and right. A series of alternating sites in between does 
not lead to an occupied cluster (case (c) of Section II) and 
at the least a two-site cluster (••) needs to be present 
(Fig. [2^1.). For the final cluster to attain length n, the 
combined lengths of the alternating sequences abutting 




FIG. 1: (Color online) Density of occupied sites: Shown are 
the density of occupied sites at the steady state, ps, (solid 
curve) and after one time step, pi, (broken line), as a function 
of the initial density, p. The initial density p is plotted as well 
(straight line), for the sake of comparison. 
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FIG. 2: Cases contributing to the eventual formation of n- 
clusters of occupied sites (see text). 



the •• seed should be 2(n — 2). Because the seed could 
be initially in any of ti — 1 locations (including right next 
to the 1=1 blocks) , the total contribution of this case is 



Pain)=q^p^{n-l)ipqr 



n> 2. 



The initial seed could be three sites long (Fig. [5]d), 
leading to a contribution 



Pb{n)=qV{n-2){pq) 



n— 3 



n > 3. 



All other cases are exemplified by Fig. [5}:. The shaded 
block near the center denotes an arbitrary sequence of 
sites with the property that it contains no two consecu- 
tive empty sites (oo). All of the sites in such a sequence 
would end up being occupied in the steady state. Denot- 
ing the probability of an m-sequence of this type by ipm, 
we can express the contribution of this last case as 

n-4 
fc=0 

To compute ipm denote by f„ the probability of an n- 
sequence that contains no oo's and ends with •, and let 
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FIG. 3: (Color online) Probability of large occupied clusters: 
The probability of n-clusters of occupied sites scales as r"^, 
in the steady state, as opposed to ~ p", of the initial con- 
figuration. Plotted are r+ (curve) and p (straight line) as a 
function of the initial density p. 



T]n denote similar sequences, but ending with o. The two 
quantities satisfy the recursion equations: 



Cn + l = {£.71 + Vn)P': Vn+1 = CnQ ■ 

Using ipn — Cn + Vm this can be recast into 



The solution is, 



c+r" + c_r" 



± 



(4) 



(5) 



where c± are determined from the boundary conditions. 

Although it is now straightforward to compute P(n) = 
Pain) + Pb{n) + Pc{n), the resulting expression is cum- 
bersome and not particularly instructive. Instead, we ob- 



serve that, for large n, P{n) 



Since r+ > p (Fig. [3]), 



the probability for large clusters in the steady state is 
exponentially larger than in the initial state. For exam- 
ple, for small p, tj^. ^ ^/p, and P{n) ~ p"/^, compared to 
~ of the initial state. 

Despite the cumbersome form of P{n), the average 
cluster size in the steady state yields a surprisingly simple 
expression; 



J2nP{n) _2-p 

j:p{n) 
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where for the last equality we are using (n,)o = 1/<Z, for 
the sake of comparison with the initial state. 

As a simple test of our results one can compute Ps = 
(n,)/((n,) -I- (?io)), where the average length of clusters 
of empty sites, (rio) = (2 — q)/p'^, is obtained from (n,) 



by exchanging p ^ q. It is easy to confirm that this 
agrees with ps of Eq. ([3]). Also, l/((n,) -I- (no)) = Kg of 
Eq. ([2]), suggesting that there is no correlation between 
the lengths of adjacent occupied and empty clusters. This 
seems remarkable, when contrasted with the correlations 
observed in P{n) for large n. 



IV. DYNAMICS 

To analyze the transient behavior of the one- 
dimensional NCO model, consider the probability ttj that 
an initially occupied site flips exactly t times before set- 
tling down for good, ttq = 2p'^—p^ is simply the probabil- 



ity that the initially occupied site is stable, i.e., it is the 
central site in one of the configurations: o»»,»»o. 
To TTi contribute the •-sites right next to the blocks i=i 
in cases (a), (b), and (c) of Section II, resulting in 
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TTi = 2p^q 



pq 



1 



q p 



. 1 +Pg 
1 



pq i^-pq 

To TT2 contribute the •'s two sites away from the blocks 
in cases (a), (b) and (d), etc. The general result is 



r, 2 2 (PlY , I 
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(7) 

pq ' -■ i~pq 

where oj = p'^q or q'^p for even and odd n, respectively. 
The analogous probabilities, at, for a o-site to flip ex- 
actly t times before stopping are obtained, as usual, by 
interchanging p -o- q in ttj . 

Armed with tt^ and at one can compute various quan- 
tities of interest. For example, the steady-state density 
is 



Ps 



E 

T=0 



{'K2t + f2r+l) 



and it is easy to confirm the agreement of this formula 
with ([3]). More importantly, the probability that a site, 
of any kind, is still active at time t (persistence), is 



St ^Y.^TTr + (Jr) - {pqY 



(8) 



The rapid exponential decay of St is typical of other dy- 
namical quantities. In practical terms, the system settles 
almost instantly to the steady state. An example of this 
effect can be seen in Fig. [U where the density of occu- 
pied sites after merely one step, pi, is compared to the 
density ps of the steady state. An exception is the proba- 
bility Pt (n) for finding n-clusters of occupied sites after t 
steps. For small p (and large n) , it converges slowly to the 
steady state, as Pt{n) ^ p"'", with at = {t+ l)/{2t+ 1). 



V. NCO ON THE BETHE LATTICE 

We now consider the NCO model on the Bethe lattice 
with coordination number z > 3 (the line can actually 
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be viewed as a Bethe lattice with z = 2). The Bethe 
lattice with z > 3 is far richer than the line, in that it 
sustains phase transitions: being infinite-dimensional, its 
phase transitions are mean-field in character, as expected 
above the critical dimension. For percolation, for exam- 
ple, a phase transition occurs at p — pc — ^/{z — 1). For 
p > Pc there exist an infinite number of infinitely large 
connected clusters of occupied sites, which are character- 
ized by mean-field critical exponents. This is known as 
the critical phase. In regular lattices there is no critical 
phase, and instead there is only one infinite cluster at 
p > Pc, known as the percolative phase Q. 

For z — 3 the neighborhood of each node consists of 4 
sites (including the site itself) and according to the NCO 
rules a site retains its own state in the event of a tie. It 
follows that, just as for the line, a two-site cluster is for- 
ever stable. At p > 1/2 the initial configuration possesses 
infinitely many infinite percolation clusters. Those clus- 
ters will remain stable indefinitely, in the NCO model, 
perhaps accreting even more of the neighboring sites as 
the updates take place. Indeed, the threshold Pc for the 
emergence of infinite clusters, in the NCO model, is ac- 
tually smaller than 1/2 (of regular percolation). 

We consider a "root" node at ^ = and define layer 
£ as the set of nodes £ links away from the root. Let 
Ni denote the number of paths, from the root to layer £, 
that never encounter an initially stable empty site, that 
is, none of the empty sites on a path have an empty 
neighbor (either on, or branching away from the path). 
All of the nodes on the Ni paths end up occupied in the 
steady state. Following a similar reasoning to the one we 
applied for ipn on the line, we find 

Ne+i^2pNi + 4p^qNi_i, 

leading to Ne ^ r^, with r = p(l -I- -y/l -|- 4g). The critical 
point for the emergence of infinite clusters in the NCO 
occurs therefore for pc « 0.344446 — the root of r = 1. 

To sum up, the NCO model on the z = 3 Bethe lattice 
is very much like percolation in that there is a critical 
phase of infinitely many infinite clusters at p > Pc, and 
only finite clusters ed p < Pc- That there is no transition 
to a percolative phase is clear from the fact that stable 
oo-clusters are seeded, initially, for all p (other than at 
the extremum, p = 1), preventing the coalescence of all 
the infinite clusters into a single one. 

The situation is even more interesting for z — A. In 
this case, the neighborhood consists of 5 sites and a min- 
imum of 2 occupied neighbors are required to stabilize 
a •-site. It is easy to see that under this condition the 
z = A Bethe lattice does not support finite stable clus- 



ters. However, the infinite percolation clusters are stable 
also in the NCO model, so it follows that the existence of 
infinitely many infinite clusters of occupied sites is guar- 
anteed for p > l/(z— 1) ~ 1/3, and may occur somewhat 
sooner, at p = pd < 1/3. At p < Pci, infinite clusters fail 
to form and the steady state is uniformly empty, since 
finite occupied clusters are unstable and vanish as the 
NCO updates evolve. The lattice becomes uniformly oc- 
cupied, at the steady state, for p > pc2 — 1 —Pci, because 
of the symmetry between empty and occupied sites under 
p ^ q. In other words, the NCO model has now three 
phases: no infinite clusters, for p < pd', infinitely many 
infinite clusters, for pd < p < Pc2', and a single infinite 
cluster (consisting of the whole lattice), for p > pc2- Un- 
like regular percolation, there are no finite clusters in the 
non-critical phases. The situation is similar for z > A. 



VI. DISCUSSION 

In conclusion, we have presented an exact solution of 
the NCO model on the line. This analysis is particularly 
interesting in view of the suspected connection between 
the NCO model and TIP Q. It should be noted that 
while invasion percolation on the line is trivially defined 
and solved, it is not at all clear how to implement the 
trapping condition — that is, superficially, there seems 
to be no difference between invasion percolation and TIP 
on the line, and neither displays the richness of results 
found for the NCO model in this paper. 

The situation is even more interesting for Bethe lat- 
tices. Just as in one dimension, there is no obvious way 
to distinguish between invasion percolation and TIP. For 
= 3, the two-phase diagram of the NCO model is the 
same as that of regular percolation (for any z). The three 
phases found for the NCO model on the Bethe lattice 
with z> A can be found in regular percolation, but only 
on hyperbolic lattices and other nonamenable one-ended 
graphs Solving the NCO model on Bethe lattices is 
an appealing open problem. It remains ultimately un- 
clear, however, how to connect the rich behavior of the 
NCO model on the line and on Bethe lattices to TIP. 
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